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Abstract 


A three-dimensional large-eddy simulation model, TASS, is used to simulate the behavior of aircraft wake 
vortices in a real atmosphere. The purpose for this study is to validate the use of TASS for simulating the 
decay and transport of wake vortices. Three simulations are performed and the results are compared with 
the observed data from the 1994-1995 Memphis field experiments. The selected cases have an atmospheric 
environment of weak turbulence and stable stratification. The model simulations are initialized with appro- 
priate meteorological conditions and a post roll-up vortex system. The behavior of wake vortices as they 
descend within the atmospheric boundary layer and interact with the ground is discussed. 


1. Introduction 

To increase airport capacity while maintaining or 
increasing the present aviation safety levels, National 
Aeronautics and Space Administration (NASA) ini- 
tiated the Aircraft Vortex Spacing System (AVOSS) 
project as an essential element of its Terminal Area 
Productivity (TAP) program to develop a real-time 
wake vortex forecasting system for Air Traffic Con- 
trol (ATC) 1,2 . 

An accurate estimation of the lifetime and trans- 
port of the wake vortices trailing behind airplanes is 
truly crucial in developing the forecasting systems 
which are responsible for aviation safety and effi- 
ciency. A focal problem is to predict the behavior 
of wake vortices under different atmospheric envi- 
ronments. Results of field experiments as well as 
numerical simulations indicate that the atmospheric 
disturbances, such as turbulence, stratification and 
wind shear have strong influence on the evolution 
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and transport of wake vortices 3 9 . 

The development of the forecasting system is 
based on both observed data and numerical simula- 
tions. During the past decades, great efforts have 
been made to understand the demise mechanism of 
wake vortices and the impact of ambient turbulence, 
stratification, wind shear and proximity effects on 
the vortices. These efforts included the laboratory 
and observed data as well as numerical modeling 9-15 . 
However, the behavior of wake vortices still remains 
elusive due to the intricate interaction of the wake 
vortex pair with its surrounding atmospheric bound- 
ary layer. On the other hand, systematic investiga- 
tion of these processes is hampered by the environ- 
mental conditions. The numerical simulations be- 
come more practical with increasing computer power 
because they allow us to study the behavior of wake 
vortices in a wider range of meteorological conditions 
than those measured during field experiments which 
are restricted to a specific location (sometimes with 
a unique climatology) and relatively short period of 
time. But numerical models must be validated using 
these observational data. 

In recent years, the large-eddy simulation (LES) 
models within the meteorological context have been 
widely used to study the behavior of wake 
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vortices 15-21 . Numerical modeling, efforts play a key 
role in the development of the AVOSS system. The 
Terminal Area Simulation System (TASS), a nu- 
merical large-eddy simulation model developed by 
Proctor 22 ’ 23 , has been adapted to accomplish this 
objective. The model has been successfully used 
to investigate the effects of atmospheric turbulence, 
stratification, wind shear and fog on the decay 
of wake vortices in both two and three dimensional 
simulations 9,17 . Proctor 9 simulated the behavior of 
wake vortices measured during the Idaho Fails and 
Memphis field programs using the TASS model. His 
results showed that the TASS model was accurate 
in simulating transport of wake vortices in two di- 
mensions for stable conditions in which there is little 
or no ambient turbulence. Han et al. 19-20 investi- 
gated the Crow instability and the effect of atmo- 
spheric turbulence on the vortex decay using three- 
dimensional TASS model. Their results agree well 
with other laboratory and field measurements. It is 
essential that all aspects of the model simulations be 
compared and validated with observed data although 
previous studies have shown that the TASS model is 
quite successful in simulating vortex decay. 

Lidar measurements, coupled with the simulta- 
neous recording of the meteorological data from the 
field measurements conducted at Memphis during 
December, 1994 and August, 1995 and at Dallas/Ft. 
Worth (DFW) Airport during 1997-1998, have pro- 
vided us with the most reliable data base so far. 
These data sets contain the information on veloc- 
ity, temperature, circulation and so on, which sup- 
ply us enough information to carefully make assess- 
ment, validation, and improvement of various numer- 
ical models. 

Ongoing studies with the TASS model are aimed 
at understanding and quantifying the decay of wake 
vortices. A determination of whether the simulations 
can realistically decay vortices is a necessary first 
step. The 'model must be validated before its results 
are used in developing predictor algorithms for an 
AVOSS system. 

The purpose of this paper is to further validate 
the TASS model. Three dimensional simulations of 
cases 1252 and 1273 from the 1995 Memphis field 
experiments are conducted. Presented are TASS re- 
sults and comparisons with the measurement data 
from both of the cases. 

In this paper, we first briefly describe the TASS 
model in section 2. We then describe the initializa- 


tion of the model and the experiments used to ini- 
tialize the model in section 3. The results from the 
simulations are presented in section 4 and, finally, 
our studies are summarized in section 5. 

2. The Model Description 

2.1 The Numerical Model 

The TASS model used in the present study is 
a three-dimensional, nonlinear, compressible, non- 
hydrostatic LES model. The TASS equation set in 
standard tensor notation is as follows 9,22 : 
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where u* is the component of velocity, t is time, 9 
is potential temperature and Oo is the reference po- 
tential temperature, p is deviation of pressure from 
the reference atmospheric pressure P, po is the ref- 
erence air density, C v and C v are the specific heats 
of air at constant pressure and volume, respectively. 
H represents the ratio of the reference density of en- 
vironment to the local density. <5*3 is the Kronecker 
delta, g is the earth’s gravitational acceleration and 
T{j and t Qj are Reynolds stress and heat flux, respec- 
tively. 


2.2 The Subgrid Turbulence Model 

The dependent variables in TASS are treated as 
averages over the grid volumes which result in sub- 
grid flux terms. The subgrid flux terms T{ 3 and tqj 
are, respectively, approximated by using first-order 
closure: 

7~ij = Po^m^dij j (3) 


T ej =p 0 K H — , (6) 
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where the deformation tensor Dij is defined as 


2 duk A 

Sxj <9x* 3 dxk tj 


(7) 


For subgrid turbulence, TASS currently uses a 
conventional Smagorinsky model with modifications 
for stratification and rotation effects: 

K m =(c s A) 2 \D\(l-aRi s -0RiR) 0 - 5 (8) 


K h = aK m , (9) 
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Here K m and Kh are the subgrid scale eddy dif- 
fusivities of momentum and heat, respectively. Ris 
is the Richardson number due to stratification, N is 
the Brunt-Vaisala frequency, A=(2Ax2Ax2Ax) 1//3 
is the filter width, and c s = 0.16 and a = 3 are con- 
stants, and (3 is a constant. The rotational Richard- 
son number, RIr, represents the rotational effect on 
turbulence and is analogous to the effect of stratifi- 
cation. 

This analogy between the rotation or streamline 
curvature and buoyancy was first suggested by Brad- 
shaw 24 . An approximate expression for RIr applica- 
ble to three-dimensional simulations is 


R*r 


jn[ 

D 2 + \D\’ 


( 12 ) 


where Tl is the magnitude of three-dimensional 
vorticity 25 . The above formulation does not require 
cylindrical coordinates or location of vortex axis, so 
it can be readily applied to fully three-dimensional 
systems. Equations (8) and (12) indicate that the 
eddy viscosity K m can be effectively reduced within 
the core of vortices when RIr or \Tt\/\D\ is large. 
Since the above formula cannot discriminate between 
a shear flow and flow with coherent rotation, a dis- 
criminator function is applied such that R{r ~ 0 when 
no coherent rotation is present (thus defaulting to the 
original Smagorinsky model for pure shear flow! 25). 
Without Eq. (12), the Smagorinsky subgrid closure 


model overpredicts the vortex decay rate in the vor- 
tex core region and causes unrealistic growth of the 
coresize 26 . 

In the TASS model, a time-splitting integration 
scheme is used, in which the time integration is bro- 
ken into small and large time steps. The acoustic 
term in the momentum and pressure equations are 
computed over the small time step, while other terms 
are calculated over the large time step. The space 
derivatives are approximated by central differences in 
quadratic conservative form with 4th-order accuracy 
for advection terms. Arakawa C grid staggered mesh 
is employed for all variables 27 . There is a choice of 
different surface boundary conditions by which one 
can specify the uniform surface temperature, the sur- 
face heat flux or the surface energy budget. There 
is also an option to choose either open or periodical 
lateral boundary conditions. A detailed description 
of the model was given by Proctor 9,22 . 

In our simulations, periodic boundary conditions 
are imposed in horizontal directions and zero gra- 
dient top boundary conditions and no slip bottom 
conditions are used in the vertical direction. Surface 
stresses are computed from similarity theory and the 
surface heat flux is prescribed. 

3. Initialization 

3.1 Vortex Initialization 


A vortex model recently developed by Proctor 28 
is used as an initial wake vortex field for three di- 
mensional simulations. This model is empirical as it 
is based on field observations of several wake vortices 
measured early in their evolution and represents the 
post roll-up, wake-vortex velocity field. The model 
represents observed wake vortex tangential velocity 
profile much better than previous models, such as 
Lamb model. Its tangential velocity, V, is repre- 
sented as 

V(r) = Ll (13) 

27cr \ / 

where Too is the circulation at r >> r c ( r c is the 
initial vortex core radius defined as radial distance 
of peak tangential velocity) and B is the span of the 
generating aircraft. The values assumed for initial 
vortex separation and circulation are derived from 
aircraft weight W, wingspan B , air density p, and 
airspeed V a , according to the conventional assump- 
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tion of elliptically-loaded wing, i.e., 


, ttB , _ 4 W 

to = — and r„ = (14) 

Note that the velocity field in Eg. (13) depends on 
the wingspan B instead of the core radius r c which 
may not be easy to measure with a good accuracy; 
Eq.(l3) is applied only at r > r c . For r < r c , the 
model is matched with the Lamb model, i.e., 

• x(l-e- 12527 ^ 2 ). (15) 

The results from Han et al . 26 show that this ini- 
tialization procedure gives a good comparison with 
observation data even though it does not depend on 
the individual aircraft configuration. 

3.2 Measurement Data 

The data we used in our simulations come from 
the experiment conducted at the Memphis Interna- 
tional Airport during August, 1995. The purpose 
of the field experiment was to measure wake vor- 
tex, meteorological and aircraft data together at an 
operational airport for use in validation of wake vor- 
tex models and in the development of the AVOSS 
prediction system. The wake vortex measurements 
were performed by a van-mounted 10.6 /zm CO 2 CW 
lidar 29 . The lidar measured line-of-sight velocities in 
a plane perpendicular to the flight path in order to 
characterize vortices generated by approaching air- 
craft. 

The atmospheric data sets were measured by sev- 
eral instruments, which included a 45m tower instru- 
mented with various meteorological sensors for wind 
velocity, temperature and humidity, a radar profiler 
and acoustic sodar for measuring winds aloft, a Radio 
acoustic sounding system (RASS) for measuring tem- 
perature data, and the radiosonde balloon launches 
for the measurement of vertical profiles of wind speed 
and temperature 30 . The standard atmospheric vari- 
ables and the atmospheric fluxes were measured at 1 
Hz and 10 Hz, respectively. 

In the Memphis field experiment, most of the tur- 
bulence data come from the tower measurements at 
5 and 40m, the turbulence kinetic energy (TKE) and 
the energy dissipation rate, e, are only available at 


heights of 5 and 40m (e was calculated by estimating 
the spectral energy density). Here, the turbulence 
intensity is defined as weak, mild or strong in terms 
of the TKE or e at a height of 40m 31 . From these 
data sets, cases 1252 and 1273 were chosen for the 
validation of the TASS model. Both cases have a 
stable stratification with weak wind shear. The val- 
ues of wind shear and stratification at the surface 
and initial vortex height for both cases are given in 
Table 1. 


Flight No. 

M-1252 

M-1273 

{dV/dz) a 

(dV/dz) h 

N\ 

N\ 

9.3 xlO ' 2 
4.0xl0 - 3 
4.1 xlO " 2 
- 6.8 xlO -4 

1.56 xlO ” 1 
-3.0xl0 ~ 3 
1 . 12 xl 0 -1 
2.50xl0- 3 


Table 1. Values of wind shear and stratification for Mem- 
phis cases 1252 and 1273 

In Table 1 , subscript “s” and “h” represent the 
values at the surface (z = 3m) and the initial vortex 
height, respectively. The vertical shear of crosswind, 
dV/dz , changes with height, and the nondimen- 

sional Brunt-Vaisala frequency is iV*=27r6o 2 /r 0 , 
where To is the initial circulation of the vortices. 

The values of the vortex parameters and flight 
data of the generating aircraft were obtained from 
the Memphis data set 29 . The flight type, initial cir- 
culation, vortex core radius, vortex separation, ini- 
tial vortex height and sink rate of vortex are shown 
in Table 2. 


Flight No. 

M-1252 

M-1273 

Aircraft 

B-757 

DC-10 

To (m 2 /s) 

345 

489 

r c (m) 

1.5 

2.5 

60 (m) 

29.8 

39.6 

h 0 (m) 

175.0 

169.5 

Vo (m/s) 

1.87 

1.97 

b o/V 0 (s) 

16.1 

20.1 


Table 2. Aircraft type and vortex parameters 


3.3 Grid/Domain Configuration 

The TASS model with the three-dimensional op- 
tion is chosen to perform all simulations in this pa- 
per. The model domains, along the axial, lateral 
and vertical directions, are 75 x 150.75 x 237 m and 
104 x 202 x 200m for the cases 1252 and 1273, re- 


4 

American Institute of Aeronautics and Astronautics 




spectively. As shown in Table 2, the initial separa- 
tion distance of the vortex pair, b 0 , is 29.8m in case 
1252 and 39.6m in case 1273. The domain size in 
the axial direction is about 2.56o for cases 1252 and 
1273. Such a small domain size in the axial direction 
may suppress the development of Crow instability 
of which theoretical maximum wavelength is about 
8.6&o 10 and cause homogeneous vortex decay along 
this direction in a statistical sense, but it can save 
much computing time without destroying the main 
vortex decay mechanism 15 . The domain sizes in lat- 
eral directions are approximately 5&o for both cases. 
The domain sizes in lateral and vertical directions 
are reasonable for a vortex pair separation of &o, be- 
cause the strain rate of neighboring periodic vorices 
is only a few percent of that exerted by one vortex 
on the other in the vortex pair. 

To save computing time, the TASS model with 
two meshes of different grid resolution is used to per- 
form our simulations. We first ran the TASS model 
with a coaxse mesh to generate a steady state tur- 
bulence. Then all variables are interpolated to fine 
mesh to simulate the vortex decay. Both coarse and 
fine meshes cover the same domain. A detailed de- 
scription of this procedure will be given in the next 
subsection. 

The grid points are, respectively, 50 x 102 x 136 in 
the coarse mesh and 50 x 201 x 316 in the fine mesh, 
with grid size (Ax, Ay , A z) = (1.5m, 1.5m, 1.5m), 
and (Ax, Ay, A z) =(1 .5m, 0.75m, 0.75m) for case 
1252. In case 1273, the grid points are 52 x 100 x 100 
and grid sizes are (Ax, Ay, A z) = (2m, 2m, 2m) 
in the coarse mesh, while the grid points are 52 x 
202 x 200 with grid sizes (Ax, Ay, A z) = (2m, lm, 
lm) in the fine mesh. The x, y and z correspond to 
the axial, lateral and vertical directions of the vortex 
system and corresponding velocity components of u, 
v and w , respectively. 

3.4 Turbulence Initialization 

To simulate the behaviour of wake vortices in a 
real atmosphere, it is of crucial importance to get an 
initial field which is analogous to real atmospheric 
environments. One of the major difficulties in initial- 
izing the model for real cases is obtaining a represen- 
tative ambient turbulence field. The difficulty arises 
from the intricate interaction of the atmospheric tur- 
bulence, stratification and wind shear. In case 1252 
and 1273, the turbulence is weak due to stable strat- 


ification, especially for case 1252. It is difficult for 
the model to generate or maintain the observed ini- 
tial turbulence field under stable stratification con- 
ditions without an external forcing because the com- 
putational domain in the model does not cover the 
whole atmospheric boundary layer and there is no 
energy source to maintain the turbulence. 

In our simulations, we first generated a turbu- 
lence field in coarse mesh through the surface heat- 
ing. After turbulence reaches a steady state, we then 
interpolated the values of all variables in the coarse 
mesh to the fine mesh. At that time, the wind and 
temperature profiles deviated from their initial val- 
ues beause of the surface heating and mixing turbu- 
lence. To maintain the initial wind and temperature 
profiles, we continued to run the model and forced 
the mean wind and temperature profiles back to their 
initial values at each time step through the surface 
cooling. During this period, we monitor the decrease 
in turbulence kinetic energy (TKE) with time at a 
height near 40m. The wake vortex pair was injected 
into the simulation when the TKE at this height was 
approximately equal to the observed TKE measured 
at 40m 

In this paper, three simulations were performed 
using the three dimensional TASS model and the 
cases 1252 and 1273 in the Memphis field measure- 
ments were used to validate the model. The first 
two simulations were designed to further simulate 
the vortex decay in case 1273 and 1252, respectively. 
The third simulation was designed to examine the 
simulation of case 1252. 

The vertical profiles of measured temperature and 
wind velocity that represent the air mass surrounding 
the wake vortices were used to initialize the simula- 
tions (Figs. 1 and 2). We assumed that the ambient 
temperature and wind velocity are homogeneous in 
horizontal directions and the ambient velocity along 
the axial direction of wake vortex (u components 
zero. 

The energy dissipation rates and turbulence in- 
tensities of the measurements and the three simula- 
tions for both case 1252 and 1273 are given in Table 
3. 

Here v =(ebo) 1/3 /V 0 and F 0 =r/27r6 0 . “ob” and 
“sm” represent the values from the observation and 
simulation, respectively. A second simulation for 
case 1252 (referred as M-1252-II) is a re-simulation 
of the first simulation, but with somewhat larger tur- 
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© y (X) u and v (m/s) ©, (K) u and v (m/s) 


Figure 1. The vertical profiles of (a) virtual potential 
temperature and (b) wind velocities for case 1252. u 
(dashed line) and v (solid line) are wind velocity in axial 
and lateral (crosswind) directions, respectively. 


flight No. 

M-1273 

M-1252 

M-1252-II 

e(°b) 
e (sm) 
V (ob) 
T) (sm) 

3.18xl0~ 4 
2.95 xl0~ 4 
1.21xl0 -1 
1.42X1CT 1 

1.35xl0 _s 
6.27x 10~ 6 
1.84xl0 -2 
3.05xlCT 2 

1.35xl0" 6 

4.30xl0~ 5 

1.84xl0 -2 

5.85xl0~ 2 


Table 3. Turbulence intensity 


bulence intensity than was observed. 

4. Results and Discussion 

4.1 Case 1273 

In case 1273 a wake vortex is generated by a DC- 
10 at an elevation of 170 m. The vortices descend 
through weak crosswind and low turbulence and are 
not influenced by the ground until after 60 s. The at- 
mospheric temperature profile is weakly stable above 
50 m elevation, but becomes increasing stable toward 
the ground. Results from this case are shown in Fig. 
3-6. 

Figure 3 shows the radial distribution of normal- 
ized circulation for 3-D simulation and Lidar mea- 
surement data for case 1273 at t= 15, 45 and 75 
second. Comparison of the profiles shows that the 
circulation computed from the model agrees gener- 
ally with that measured by the Lidar. Both the sim- 
ulations and observations show that the circulation 
decreases with time and the core radius remains a 
constant. 


Figure 2. Same as Fig. 1, but for case 1273. 

The circulation decays more or less uniformly with 
radius as shown in Fig. 3. At t=75s, the vortex has 
descended to about 35m above the ground (see Fig. 
6), where the stratification becomes more stable and 
the ground begins to interact strongly with the vor- 
tex pair. 

Figure 3 also indicates that the modeled circu- 
lation, in some degree, deviates from the observed 
circulation. One of the reasons is that the mod- 
eled circulation was obtained directly by integrating 
the vorticity field and then averaging the circulations 
at all cross-sections along axial direction, while the 
measured circulation was estimated from Lidar line- 
of-sight velocities at one cross-section. 

The radial distribution of tangential velocity from 
the simulation and Lidar measurement at t=15, 45 
and 75s is shown in Figure 4. It is clear that the tan- 
gential velocity decays with increasing time and ra- 
dial distance. At t=75s, the tangential velocity near 
the edge of the core radius is reduced to about 30% 
of its initial value. The rapid decay is caused by the 
ambient turbulence and stratification, rather than by 
molecular diffusion and core expansion which need 
more than several hours to undergo a similar decay. 
In fact, Fig. 4c shows that the core radius only in- 
crease to 3.0m from its initial value of 2.5m. These 
results agree well with the estimates of Sarpkaya 31 . 

The averaged circulations of 3- 10m and 5- 15m are 
plotted in Figure 5. The 3-10m averaged circulation 
is well matched with measurements, except for the 
time less than 10s when the observed circulation was 
larger than the modeled circulation. In our simula- 
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Circulation Profiles for 3-0 Case 1273 ( t = 15 sec ) 



Circulation Profiles for 3-0 Case 1273 ( t = 45 sec ) 



Circulation Profiles for 3-0 Case 1 273 ( t = 75 sec ) 



Figure 3. The radial distribution of circulation from the TASS simulation and Lidar measurements for case 1273 at 
(a) t=15s, (b) 45s and (c) 75s. Circulation is normalized by the initial theoretical circulation, and radius is normalized 
by the initial core radius (rc=2.5 m). Symbols represent Lidar measurements of left and right sides of port vortex. 
Solid and dashed lines represent the simulated Port and Starboard vortices, respectively. 
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Circulation Profiles for 3-0 Case 1273 ( t = 15 sec ) 


Circulation Profiles ior3-D Case 1273 ( t » 45 sec ) 




Circulation Profiles tor 3-0 Case 1273 ( t = 75 sec ) 



Figure 4. The radial distribution of tangential velocity from the TASS simulation and Lidar measurement for case 
1273 at (a) t=15s, (b) 45s and (c) 75s. Symbols represent Lidar measurements of left and right sides of port vortex. 
Solid and dashed lines represent the simulated port and starboard vortices, respectively. 
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t (s) 1 (s) 


Figure 5. The 3-10m and 5-15m average circulations versus time from both the simulation and measurements for 
case 1273. (a) the 3-10m average circulation; (b) the 5-15m average circulation. 



Figure 6. The vortex height history from the simulation and field measurements for case 1273. Symbols represent 
Lidar measurements. Solid and dashed lines represent the simulated port and starboard vortices, respectively. 
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tion, the 5-10m averaged circulation decays slowly 
before t=60s and more rapidly after that time (Fig. 
5b). Although the modeled circulations are slightly 
larger than the observed circulation, the differences 
are within 20%. The discrepancy at early times may 
be caused by measurement error, since Campbell et 
al. 32 have found that the measured Lidar circulations 
were overestimated during initial times. 

Comparisons of the modeled results with field 
measurements of the vortex trajectories are shown 
in Figure 6. It can be seen that the development 
of the vortex height with time from both the simu- 
lation and observation is in reasonably good agree- 
ment. The vortices linearly descend until at t=60s 
when they stall due to the stable stratification, wind- 
shear, and ground effect. The vortices descended to 
an average height of 55m at about t=85s. After that 
time, they rise slightly. 

As shown in Fig. 6, the port vortex is deflected 
to a higher altitude than its starboard counterpart, 
and the mutual induction between the two vortices 
is weakened because of the increasing separation be- 
tween the two vortices. Our simulation is consistent 
with Proctor’s two-dimensional simulation 17 , regard- 
ing the influence of crosswind shear on vortex de- 
scent. Compared with the measurement data (Fig. 
6), the starboard vortex descended a lower altitude 
(about 25m) than that observed (about 40m) during 
later time of the experiment. 

4.2 Case 1252 

Case 1252 is somewhat similar to 1273 except 
that the wake is generated by a B-757 in an environ- 
ment of slightly weaker stratification and turbulence. 
Both the simulated and observed wakes descend to 
the top of the stable layer at 40m in about 100s. 

The radial distribution of normalized circulation 
and tangential velocity from the 3-D simulation and 
Lidar measurement data for case 1252 at t= 15, 30 
and 80s are shown in Figures 7-8. The modeled circu- 
lation decays slower than the measured circulation, 
particularly at t<40s or r/rc > 5. The tangential 
velocity is well matched with the observation data, 
except near the core radius where the model under 
predicts the tangential velocity. 

Figure 9 shows the averaged circulations of 3- 10m 
and 5- 15m at different times. The model predicts the 
decay of the 3-10m averaged circulation very well, 


but under predicts the decay of the 5- 15m averaged 
circulation. Measurement shows decay more rapidly 
before t=40s. One possible reason for the under pre- 
diction of the decay could be an underestimation in 
the intensity of the ambient turbulence. Therefore, 
an additional simulation has been run using a slightly 
higher value of turbulence and will be presented later 
in this section. 

The vortex trajectories from the simulation and 
observation are shown in Figure 10, and it can be 
seen that the vortex trajectories are well predicted. 
Unlike in case 1273, the vortices do not show any 
rising motion prior to ground interaction. Fig. 10a 
shows the wake vortices linearly descending with a 
velocity of 1.8m/s before 70s, but slowly descending 
after 70s apparently due to the strong stratification 
and ground effects. The effect of the crosswind is 
clear in Fig. 10b. Both vortices are transported in 
the same direction. They move away from the posi- 
tion of generation before 45s and turn back to their 
original position during later time. This is partly 
due to the downward transport of crossflow momen- 
tum from their generation point. It is worthy to note 
that the effect of the crosswind is not simple advec- 
tion since the slope of the lateral position vs time 
relation is different for both vortices. 

To investigate whether the relatively faster decay 
from the Lidar measurements might have been as- 
sociated with an under-reported level of turbulence 
intensity, we have rerun case 1252 with a slightly 
larger value for ambient turbulence. Otherwise, all 
conditions are identical with those in the previous 
simulation, with e now set to 4.3 x 10“ 5 m 2 /s 3 . 

The radial distribution of normalized circulation 
and tangential velocity for the simulation and Li- 
dar measurement data at t= 30 are shown in Figure 
12. It can be seen that the modeled circulation and 
tangential velocity are well matched with the Lidar 
measurements. Comparing with Figs. 8-9, the circu- 
lation decay is much faster at all radii for the larger 
intensity of ambient turbulence. 

Figure 12 shows that the averaged circulations of 
both 3-10m and 5-10m are well predicted except that 
the model slightly over-predicts the decay of the 3- 
10m averaged circulation after t=60s due to slightly 
stronger turbulence. Compared with Fig. 9, we can 
see from Fig. 12 that this simulation better repro- 
duced the behavior of wake vortices. The vortices 
decay rapidly for the time less than 40s and decay 
slowly after that time. A possible explanation for 
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1.1 


Circulation Profiles for 3-0 Case 1252 { t = 1 5 sec ) 



Figure 7. The radial distribution of circulation from the TASS simulation and Lidar measurements for case 1252 at 
(a) t=15s, (b) 30s and (c) 80s. Circulation is normalized by the initial theoretical circulation, and radius is normalized 
by initial core radius (rc=1.5 m). Symbols represent Lidar measurements of left and right sides of port vortex. Solid 
and dashed lines represent the simulated port and starboard vortices, respectively. 
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Figure 9. The 3- 10m and 5- 15m average circulations versus time from both the simulation and measurements for 
case 1252. (a) the 3-10m average circulation; (b) the 5-15m average circulation. 
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Figure 10. The vortex trajectories from the simulation and Lidar measurements for case 1252. (a) vortex height 
history; (b)vortex lateral position. Symbols represent Lidar measurements. Solid and dashed lines represent the 
simulated port and starboard vortices, respectively. 


Circulation Profits tor 3-0 Casa 1252 ( t « 30 sec ) 



Tangential Velocity Profiles for 3-D Case 1252 { t = 30 sec ) 



Figure 11. The radial distribution of (a) circulation and (b) tangential velocity from the TASS simulation and Lidar 
measurements at t= 30s. Circulation is normalized by the initial theoretical circulation, and r c is initial core radius. 
Symbols represent Lidar measurements of left and right sides of port vortex. Solid line represents the simulated port 
vortex. 
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Average Circulation Evolution for 3-0 Case 12S2 ( r » 5 — 15 m ) 


Aveng* Circulation Evolution ter 3-0 Case 1252 ( r - 3—10 m ) 




Figure 12. The 3- 10m and 5- 10m average circulations versus time from the simulation and Lidar measurements, (a) 
the 3- 10m average circulation; (b) the 5- 10m circulation. 


this good agreement is that the measured dissipation 
rate was under-estimated or that the value observed 
at 40 m did not characterize the turbulence intensity 
at other altitudes. 

5. Summary 

To validate the application of LES to wake vor- 
tex decay, three-dimensional simulations in a realis- 
tic atmospheric boundary layer were performed us- 
ing NASA’s TASS model. Two case studies were 
presented which were initialized from and compared 
with Memphis fields measurements. 

Our results indicate that the TASS model suc- 
cessfully predicted the behavior of aircraft wake vor- 
tices within environments representing a transition- 
ing atmospheric boundary layer following sunset. 
Good agreement was obtained between simulations 
and measurements for the wake vortex trajectories. 
Exceptional agreement with the measurements was 
obtained for vortex decay. The three-dimensional 
simulations suggest that ambient turbulence, stratifi- 
cation and wind shear play important roles in the de- 
cay and transport of wake vortices. Ambient turbu- 
lence significantly influences the vortex decay. Study 
on the specific roles of wind shear, stratification and 
turbulence on vortex behavior pertinent to these two 
cases are in progress. 
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